Recent results on D° mixing from Belle 
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We report on recent measurements of the D° — D° mixing and CP violation parameters performed 
by the Belle experiment. The evidence for the mixing phenomena in the system of neutral D mesons, 
arising in the study of D° — > K + K~ ,tt + tt~ decays is presented first. Using a time dependent Dalitz 
analysis of D° — > A"g7r + 7r~ decays we also obtained the most precise up-to-date determination of the 
mass difference of the two D meson mass eigenstates. The presented results are based on 540 fb _1 of 
data recorded by the Belle detector at the KEKB e + e~ collider. We conclude with short prospects 
for the future measurements. 



INTRODUCTION 

Measurements in the field of charmed 
hadrons experience a revival in the recent 
years. The reason for an increased inter- 
est is twofold: the B-factories provide for an 
abundant source of charmed hadrons. The 
integrated luminosity L « 700 fb _1 of the 
KEKB collider [ij corresponds to a produc- 
tion of around 900 x 10 6 charmed hadron 
pairs in a clean environment of e + e~ col- 
lisions. Secondly, a dual role of charm 
physics is exploited: as an experimental 
test ground for different theoretical predic- 
tions, most notably the lattice QCD, en- 
abling in turn a more precise determinations 
of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements; and as a standalone field 
of various Standard model (SM) tests and 
searches of new physics (NP) phenomena. 

Search for the D° — D mixing (a quest 
started soon after the discovery of Z? mesons 
in 1976 01) belongs to the latter category. 
It is governed by the lifetime of D mesons, 
t = 1/r, and by the mixing parameters x = 
(m 1 - m 2 )/r and y = (Fi - T 2 )/2T. m 12 
and r 12 denote the masses and widths of the 
mass eigenstates D\ and D 2 , respectively, 

\D h2 )=p\D }± q \lf) . (1) 

T = (Ti + T 2 )/2 is the average decay width. 
The mixing rate is severely suppressed due 



to the small SU(3) flavor symmetry break- 
ing (ml ss to 2 d ) and smallness of the \V u b\ 
CKM matrix element Q . Calculations based 
on the effective AC = 2 Hamiltonian (i.e. 
contribution of the box diagram, providing 
a satisfactory description of mixing in the 
systems of K°, B® and B® mesons) yield a 
negligible mixing parameter magnitude \x\ ~ 
0(1O~ 5 ). Long distance contributions to the 
D - D° transitions are difficult to calculate. 
Approaches based on the operator product 
expansion Q or summation over exclusive 
intermediate states, accessible to both D° 
and D° 0], result in \x\,\y\ < 10 -3 and 
Mi \y\ < 10~ 2 , respectively. 
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FIG. 1: Relation between the mixing parameter 
x and product of CKM matrix elements in model 
including the fourth generation down-like quark 
b' @. 

Regardless of large uncertainties of the 
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predictions, the SM probability of D° meson 
to oscillate into its antiparticle before decay- 
ing, R M w (x 2 + y 2 )/2, is small, at most of 
the order 0(1O -4 ). New as- yet-unobserved 
particles could contribute to the loop di- 
agrams and thus significantly increase the 
value of \x\. A limit of \x\ < 0(1O~ 2 ) would 
put significant constraints on the parameter 
space of large number of NP models @ . As 
an example a possibility of the fourth family 
of fcrmions, including a down-like b' quark 
of mass m&', may be considered (Fig. [!}. 
An experimental value of \x\ < 10 -2 would 
constrain the CKM matrix elements product 
\V uV V cb ,\ < 3 x 1(T 3 for m b , > 200 GeV/c 2 , 
a limit much more stringent than the one 
following from the current CKM matrix uni- 
tarity requirement. Similar parameter con- 
straints can be put to 17 out of 21 NP models 
considered in [6(. The width difference (y), 
on the other hand, is governed by D decays 
into physical states, where no significant de- 
viations from the SM have been observed up 
to date. However, y is known to vanish in the 
limit of exact SU(3) flavor symmetry. Hence 
NP contributions not vanishing in this limit 
could affect the value of y regardless of their 
small contribution to the decay amplitudes 



Beside the possible effects on the mixing 
parameters, NP could produce a sizable vi- 
olation of the CP symmetry (CPV) in D 
meson decays 0. Within the SM the CPV 
is expected to be small. Of the three types of 
the CP violation, CPV in decays, CPV in 
mixing and in the interference between mix- 
ing and decays, the first one is expected to 
be present only in the singly Cabibbo sup- 
pressed decays. To these, beside the tree 
amplitude also penguin diagrams (c — > uqq) 
can contribute, and the existence of at least 
two amplitudes of different strong and weak 
phase is a necessary condition for this type 
of violation to occur [s|. The weak phase 
difference between the two amplitudes is ~ 
^{VcdV* d V cs V* s ) < 10" 3 , which represents 
a rough estimate of the expected CPV ef- 
fect [3J. The asymmetries due to CPV in 
mixing and interference can be expressed as 



Acp cx ycoscp and oc xsincp, respectively, 
where <f> ~ 10~ 3 is the weak phase between 
the mixing and decay amplitudes. Hence the 
magnitude of these types of CPV is even 
smaller. Observation of the CP violation an 
order of magnitude larger than these expec- 
tation would clearly point to the intervention 
of NP. 

The presented studies were carried out 
by the Belle detector, a general purpose 
full solid-angle spectrometer (T3] operating 
at the asymmetric e + e~ KEKB collider 1]. 
The center-of-mass (CMS) energy of the col- 
lisions corresponds to the mass of the T(4S), 
decaying to a pair of B mesons. In addi- 
tion to the BB production, the cross sec- 
tion for the continuum production of u, d, s 
and c quark pairs through a virtual pho- 
ton exchange at this energy is several times 
larger. In the presented measurements the 
D mesons produced in e + e~ — * cc are 
reconstructed [2^]. Two main features of 
the detector are exploited for this purpose. 
The identification of detected charged tracks 
is performed using a combined information 
from several detector sub-modules 111. An 



illustration of the performance can be given 
by the efficiency for the charged kaon iden- 
tification, ~ 90%, with the tt^ 1 misidentifica- 
tion rate < 10% for tracks with momenta be- 
tween 1 GeV/c and 3.5 GeV/c. A silicon ver- 
tex detector enables a precise determination 
of the decay time t of short-lived particles. 
For D° — > K~tt + decays, the distribution of 
estimated uncertainties on t peaks at around 
t(D°)/3 and has an average of t(D°)/2 (see 
Fig. left 0); t(D°) is the world average 
value of the D° lifetime Ftl. 



MEASUREMENTS 

There is a long list of measurements de- 
voted to the D°- D° mixing from the Belle 
collaboration. Both, semilcptonic [3] and 
hadronic decays [IH have been exploited in 
the past. The most sensitive recent mea- 
surements, using decays to a CP eigenstate 
(D° -> f C p, with f C P = K+K-, tt+tt-) [13 
and to a self-conjugate final state (D° — > 
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FIG. 2: Left: Normalized distribution of esti- 
mated errors on the decay time t in D° — > K~tt + 
decays. Fraction fa of events has an uncer- 
tainty of (Tj. Right: Measured D° lifetime from 
D° — > K~% + in different running periods. The 
value shown on the top is the average. The left- 
most point shows the current world average of 
lifetime fTl. 



KgTT + ir ) flq . are presented in this paper. 

There are several methods and selection 
criteria in common to the presented measure- 
ments. In order to search for events where a 
D° undergoes a transition to D the flavor of 
the initially produced neutral meson must be 
tagged. This is achieved by reconstruction of 
decays D*+ -> P°tt+ or D*~ -> ifn-. The 
charge of the characteristic low momentum 
pion ir s tags the flavor of the initially pro- 
duced D meson. The energy released in the 
D* decay, 



M(D*) - M(D°) - m„ 



(2) 



has a narrow peak for the signal events 
and thus helps in rejecting the combinato- 
rial background. Here, M(X) is used to 
denote the invariant mass of the X decay 
products, and mx for the nominal mass of 
X. D° mesons produced in B decays have 
different decay time distribution and kine- 
matic properties than the mesons produced 
in continuum. In order to obtain a sample of 
neutral mesons with uniform properties we 
require the momentum of the reconstructed 
D* mesons in the CMS to be larger than 
2.5 GeV/c 2 . Since the momentum of D* in 
B — > D*X is kinematically constrained, this 
requirement completely rejects the D 0, s from 
the latter source. Last but not least, the se- 
lection criteria are optimized using the MC 
simulation, in order not to bias the results of 
the measurements. 



Evidence for charm mixing in 

D° -> K+K-/TT+1T- 



In the limit of no CPV the mass eigen- 
states of neutral charmed mesons, with dis- 
tinct values of lifetime l/ri.2, are also CP 
eigenstates. Hence only the mass eigen- 
state component of D° with the CP eigen- 
value equal to the one of fcp contributes 
to D° — > fcp decays. By measuring the 
lifetime of D° in decays to fcp one deter- 
mines the corresponding 1/ri or I/T2. On 
the other hand, both CP states contribute 
in decays to non-CP final states, like K~tt + . 
The measured value of the effective lifetime 
in the latter process corresponds to a mix- 
ture of 1/Ti and l/IV By explicit writing 
of decay time rates in the presence of oscilla- 
tions, and taking into account \y\ << 1, one 
derives a relation between lifetimes as mea- 
sured in D° — > fcp and in decays to a mixed 
CP final state to be [ltI 



r(fcp) 



r(D°) 



1 



VfVcp 



(3) 



with rjf = ±1 denoting the CP eigenvalue 
of fcp- T (D°) represents the effective D° 
lifetime as measured in decays to non-CP 
eigenstates (1/r) and the relative difference 
of the lifetimes is described by the parameter 
Vcp- 

The final states fcp = K + K~ ,it + tt~ arc 
CP eigenstates with T)k+k- tt+tt- = +1- 
The ratio of lifetimes measured in these de- 
cays and in the D° — > K~ir + yields the value 
of yep- 



t(K ir + ) 
Vcp = > - 1 



(4) 



r(fcp) 

Expressed in terms of the mixing parameters 
yep reads [T3| 

1 



Vcp = V cos 4> - -A M sin < 



(5) 



with Am and 4> describing the CPV in mix- 
ing and interference between mixing and de- 
cays, respectively. If for the moment the 
possibility of CPV is neglected (Am, 4> = 0; 
search for the CPV is described separately 
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TABLE I: Signal yields and purities of selected 
samples. 



Final state 


Signal yield 


Purity 


K+K~ 


lllxlO 3 


98% 




1220x10^ 


99% 


7T + 7T~ 


49xl0 3 


92% 



in a later section), one notes that yep = V- 
The described method of yep determination 
has been exploited in 

Selection of the D° candidate decays is 
based on M(D°)/<tm (where om is the de- 
cay channel dependent resolution), Aq = 
q — (mi). — mo — Tn„) and at- Distributions 
of M(D°) and q for D° — > K + K~ are shown 
in Fig. [3J The signal yields and purities of 




1.82 1.84 1.86 1.88 1.9 0.005 0.01 0.015 0.02 
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FIG. 3: Left: M(D°) distribution for D° -» 
K+K~ with |Ag| < 0.80 MeV. Right: Aq distri- 
bution for decays with M(D°)/a M < 2.3. Full 
histograms represent the results of the fits. 

selected fcp and K~tt + samples, following 
from the fit to the tuned simulated samples, 
are given in Table [J The amount of back- 
ground in all three reconstructed channels is 
low. 

Final state tracks are refitted to a com- 
mon D decay vertex. The production point 
is found by constraining the D° momentum 
vector and tt s from the D* decay to origi- 
nate from the e + e~ interaction region. The 
proper decay time t is calculated as a projec- 
tion of the vector L joining the two vertices 
onto the momentum of D , t = m^L ■ p/p 2 . 
To determine yep we perform a simultane- 
ous binned likelihood fit to the decay time 
distributions in the three decay modes, with 
lifetimes related by a free parameter yep- 

The t distributions are described as a sum 
of the signal and background contribution 



B(t). The signal contribution is a convolu- 
tion of an exponential and a detector resolu- 
tion function R(t): 

N ■ f 

dN/dt = -5£ / e -t /r • R(t - t') dt' + B{t). 

(6) 

The composition of the resolution function 
is illustrated in Fig. [2j left. Normalized 
distributions of estimated at, based on the 
uncertainties of the decay length determina- 
tion, are plotted for individual decay chan- 
nels. In an ideal case each <Tj value represents 
a Gaussian resolution term with a weight 
/j. Study of the normalized residual distri- 
butions, (^reconstructed — ^generated )/""*) how- 
ever, reveals that they cannot be described 
by a single Gaussian function. They are well 
described by the sum of three Gaussians, 

— 1 (^reconstructed ^generated! t^fc to), 

with weights Wk , widths ct£ u11 and a common 
mean to. It follows that each Oi represents 
a resolution term composed of three Gaus- 
sians. The final parametrization of the reso- 
lution function is thus 

71 3 

R{t-t') = Y^fi^w k G{t-t';(Ji k ,t ) , 

i=l k=l 

(7) 

with (Tjj. = Sfccr^ u <7j. The scale factors 
are introduced to describe small differences 
between the simulated and real cr^ u11 . 

The background distribution is repre- 
sented by a sum of an exponential and S func- 
tion, convolved with the resolution function 
parametrized as above. Parameters of B{t) 
are determined from fits to t distributions of 
events in the M(D°) sidebands. 

Several running periods, coinciding with 
changes to the detector, were identified based 
on the resulting t(K~it + ). For one of the pe- 
riods the resolution function ([7]) is modified 
to be slightly asymmetric in order to yield a 
consistent lifetime. This is achieved by intro- 
ducing a decay mode dependent difference of 
to's of the first two Gaussian terms in R(t). 
This behaviour has been reproduced by gen- 
erating a special MC sample which includes 
a small additional misalignment between the 
vertex detector and central drift chamber of 
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the Belle detector. The lifetime measured in 
D — ► K~ir + decays shows a good consis- 
tency among the running periods as well as 
with the world average value [l3[ (Fig- El 
right). 

Simultaneous fits to decay time distribu- 
tions of K + K~, K~ir + and tt + it~ were per- 
formed for individual running periods and 
the resulting yep values averaged to obtain 
the final result. Fits are presented in Fig. 
HJa)-(c) by summing the data points and 
the fit function values. The agreement of 
the fit function with the data is excellent, 
X 2 /n.d.f = 312/289. The same is true for 
all individual fits as well. The final value 
obtained is 



t distribution of reconstructed MC events) 
and variation of selection criteria (effect es- 
timated using high statistics MC samples). 

The resulting yep deviates from the null 
value by more than 3 standard deviations 
(more than 4 standard deviations consider- 
ing the stat. error only) and represents a 
clear evidence of D° — D mixing, regardless 
of possible CPV. The difference of lifetimes 
is made visually observable by ploting the ra- 
tio of decay time distributions for decays to 
f C p and K~n + in Fig. BJd). 

Measurement of charm mixing 
parameters in D° — + Ksit + TT~ 



yep = (1.31 ± 0.32(stat.) ± 0.25(syst.))% . 

(8) 




-2000 2000 4000 2000 4000 



t (fs) t (fs) 

FIG. 4: (a)-(c): Result of the simultaneous fit 
to decay time distributions in D° decays to indi- 
vidual final states. The hatched areas represent 
the contribution of backgrounds, (d): Ratio of 
D° — > fcp and D° — > K~tv + decay time dis- 
tributions. The slope visualize the difference of 
effective lifetimes. 

The largest systematic uncertainties fol- 
low from the assumption of equal to for dif- 
ferent decay channels (estimated by relaxing 
this constraint), possible deviations of accep- 
tance dependence on decay time from a con- 
stant (estimated by a fit to the generated 



To a hadronic multi-body final state sev- 
eral intermediate resonances can contribute. 
In a specific example of the self-conjugated 
mode D° — ► A's7r + 7r~ contributions from 
Cabibbo favored decays (e.g. D° — > 
K*~ir + ), doubly Cabibbo suppressed decays 
(e.g. D° -> K*+ir-) and decays to CP 
cigenstates (e.g. D Q —> p°Ks) are present. 
Individual contributions can be identified by 
analyzing the Dalitz distribution of the de- 
cay (see Fig. [5]). 

Different types of intermediate states ex- 
hibit also a specific time evolution (due to 
their specific superposition from the D mass 
cigenstates). While the decay time distri- 
bution of CP eigenstates depends on the pa- 
rameter y 7 the t evolution of doubly Cabibbo 
suppressed decays depends on x' — x cos 5 + 
ysirtS and y' = ycosS — xsmS, where 8 is 
a strong phase difference between these and 
the corresponding Cabibbo favored decays. 
Since in decays to Kstt + it~ both types in- 
terfere it is possible to disentangle the rela- 
tive phase by performing a fit to the Dalitz 
distribution. This in turn enables a direct 
determination of the mixing parameters x 
and y instead of their rotated versions x' and 
y' . The method was successfully exploited in 

The two Dalitz variables are defined as 
m 2 _ = M 2 {K s tt-) and m\ = M 2 (K s tt + ). 
Decay time dependent matrix element for 
D° — > A's7r + 7r _ decay, where the initially 



6 




--S 3__„ 

m?(Ge\T*n 



FIG. 5: Illustration of the decay time dependent 
Dalitz analysis. Contributions of different inter- 
mediate states (sketched by white lines) can be 
identified by analysis of the Dalitz distribution. 
Different types of decays exhibit a different de- 
cay time propagation. By studying the time evo- 
lution of the Dalitz distribution one effectively 
determines the t distribution of different types 
of decays and by this the mixing parameters x 
and y. 

produced D meson is a D°, is written as 

M(m 2 _,m% t t) = (K s n + TT-\D a (t)) = 
= }-A{m 2 _,m\)[e-^ t + e-^ t ] + 

+ ^(m 2 _,m 2 + )[e^ Alt ~ e - 4A2t ] 

(9) 

In the above expression A(m 2 _, m+) and 
A(m_ , m 2 ^) are the instantaneous ampli- 
tudes for D° and D° decays. The depen- 
dence on the mixing parameters arises upon 
squaring the matrix element in which A1.2 = 
mi,2- iTi^/2. M{rn 2 _,rn 2 +1 t) describing the 

decay of an initially produced 15° is written 
in an analogous form. Neglecting CPV one 
finds M(m 2 _,m 2 ^, t) = M.(m\,m 2 _,t). 

Amplitudes for D decays are parametrized 
in the isobar model as a sum of Breit- 
Wigner resonances and a constant non- 
resonant term: 

A{m 2 _,m%) = ^a r e^"B r (m 2 _,m 2 F )+a N Re^ 

r 

(10) 



Functions B r are products of Blatt- 
Wcisskopf form factors and rclativistic 
Breit-Wigners (H)]. The described signal 
distribution is convolved by the detector 
mass resolution function (for 7r + 7r~ invariant 
mass only) and multiplied by (m_ , m 2 + ) 
dependent efficiency. The expected de- 
cay time distribution is convolved with a 
resolution function described by a sum of 
three Gaussians with a common mean. The 
mean and scale factors for the widths of the 
resolution function are free parameters of 
the fit. 




FIG. 6: Top left: 95% C.L. region for parameters 
x and y as obtained in D° — > Kstt + tv~ decays. 
Top right and bottom: Projections of the Dalitz 
distribution to squares of two-particle invariant 
masses. Full line is the result of the fit. Note 
that in the plots m± corresponds to M 2 {Ks'k ± ) 
for D° decays and to M 2 (Ksir^) for D decays. 

The isolation of signal is based on M(D°) 
and q variables, described in the previous 
section. Selected sample of decays used for 
the measurement consists of 534 x 10 3 sig- 
nal decays with a purity of 95%. Fractions 
of individual backgrounds are obtained from 
the two-dimensional fit of M(D°) and q dis- 
tributions. 

The Dalitz distribution of combinatorial 
r background (4%) is obtained from events in 
the M(D°) sidebands. This probability den- 
sity function (p.d.f.) is multiplied by a sum 
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of an exponential and 6 functions to repre- 
sent the t distribution, and convolved with 
the same resolution function as for the sig- 
nal. The background with a true D° and a 
random slow pion represents 1% of the sam- 
ple. The p.d.f. of this background is the 
same as the one of the signal. 

An unbinncd likelihood fit is performed to 
distribution of events in the signal region. 
Results of the fit in which we neglect pos- 
sible CPV are projected to the Dalitz vari- 
ables in Fig. [6l The Dalitz model which 
includes 18 intermediate states represents a 
good description of the data. 

The effective value of the D° lifetime fol- 
lowing from the fit, t(D°) = (409.9 ±0.9) fs, 
is in good agreement with the world average 
value The decay time projection of the 
fit is presented in Fig. [7] 

io \ i 



into two categories: uncertainties related to 
the Dalitz model, and others. The former 
are estimated by repeating the fit with the 
K-matrix parametrization of the scalar res- 
onances, and by estimating possible biases 
in the ratios of doubly Cabibbo suppressed 
and Cabibbo favored decays using the sim- 
ulation. The largest uncertainty in the lat- 
ter category arises from the variation of the 
selection on the D* CMS momentum and 
from the assumption of factorization of the 
Dalitz and decay time distributions of com- 
binatorial background (estimated by using 
Dalitz distribution of combinatorial back- 
ground from different t intervals). 

The resulting value of the mass difference 
x is the most precise measurement of this pa- 
rameter, improving the precision of the pre- 
vious measurement [l9j by an order of mag- 
nitude. 
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FIG. 7: Top: Decay time distribution of selected 
D° -> A'stt+tT . Full line is the result of the fit. 
Contribution of background is presented by the 
lower line. Time scale is given in fs. Bottom: 
Residuals between the data and the fitting func- 
tion. 



Finally, the result for mixing parameters 



is 



x = (0.80 ± 0.29(stat.) + °' 13 (syst.))% 
—0.16 



(11) 



y = (0.33 ± 0.24(stat.) + °' j° (syst.))% 

(12) 

The systematic uncertainties are divided 



Search for CPV 

The value of yep measured in D — ► fcp 
depends on the parameters describing CP vi- 
olation (eq. (Op). The definition of the pa- 
rameters is II 711 



IC = i 

\p\ 2 - 

qA(D° -> K+K') 
pA(D° — ► K+K ) 



Mi* 
IpI 



(13) 



where Am 7^ is a sign of the CPV in mix- 
ing and (j) ^ of the CPV in the interference 
between mixing and decay. In D decays to 
a, CP eigenstate one can define a, CP asym- 
metry: 



Ar = 



r(D° 



fcp) - r(D° -» fcp) 



r(D° 



fcp) + r(D° -> fcp) 



, (14) 

which can be expressed in terms of funda- 
mental parameters as 



Ar = ^A M y cos< 



x sin ( 



(15) 
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Hence by separately measuring the lifetime 
of D° and D tagged decays, we measure [l2[ 

A T = (0.01 ± 0.30(stat.) ± 0.15(syst.))% . 

(16) 

Systematic uncertainty receives similar con- 
tributions as in the case of yep measure- 
ment. The CPV in mixing and interfer- 
ence is thus not observed with a sensitivity 
of - 0.35%. 

To search for CPV in Ksir + n~ decays a 
more general fit than the one described in 
the previous section is performed flil ]. In 
addition to previous parameters we allow for 
\q/p\ 7^ 1 and 0^0. To check for a possibil- 
ity of CPV in decays the parameters a r and 
4> r of Eq. (fT0|) are allowed to be different for 
D° and D decays. The resulting Dalitz pa- 
rameters arc consistent for the two samples 
and no sign of CPV in decays is observed. 
Results of the consequent fit assuming no di- 
rect CPV are 



M =0 .86- 
\P\ 



-0.30 
-0.29 



(stat.) 



-0.10 
-0.09 



(syst.) 



i_n oft 

4> = (-0.24 ' (stat.) ± 0.09(syst.)) rad . 

U .o -L 

(17) 

Also this measurement shows no evidence of 
CPV. 

The 95% confidence level region in x, y 
plane, following from the no-CPV and CPV 
allowed fits to D° — > A'stt+tt- decays, is 
shown in Fig. [6j top left. It should be noted 
that if the CPV parameters are left free 
in the fit, the solution (—x,—y,arg(q/p) + 
7r) is an equally probable solution as the 
{x,y,arg(q/p)). In terms of Fig. [6] this 
means that contours reflected over the (0,0) 
point also represent an allowed region of pa- 
rameters space. A peculiar shape of the 
CPV allowed contour in the vicinity of the 
(0, 0) point is a consequence of no possible 
CPV in mixing or in interference between 
mixing and decays when x, y = (in rest 
of the region the sensitivity of measurement 
is spread among two mixing and two CPV 
parameters; close to x, y = the sensitiv- 
ity of the measurement is mainly to x and 



y and thus the likelihood function becomes 
steeper) . 

Belle also obtained a preliminary result of 
t- and Dalitz plane-integrated CPV search in 
D meson decays to tt+7t-tt°. The amount of 
signal and backgrounds in the selected sam- 
ple is determined by a fit to M(D°) and 
AI(D ) distributions, shown in Fig. [8] 
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FIG. 8: Fits to M(D°)(M(D a )) distributions of 
D°(D°) -> tt+tt-tt decays. 

The number of D decays is calculated from 
the Dalitz distribution of events by subtrac- 
tion of the background. Background distri- 
bution is determined using MC simulation 
and the uncertainty due to the modelling is 
included in the systematic error. The back- 
ground subtracted yield is corrected for the 
efficiency in bins of Dalitz plane. By per- 
forming the described calculation separately 
for D° and D tagged decays, we obtain 



A CP = 



r(D° -> TT+TT-TT ) - T(D 



T(D° -> TT+TT-TT ) +T(D 



= (0.43 ± 0.41(stat.) ± 1.31(syst.))% 



(18) 



OUTLOOK AND SUMMARY 

Belle has recently presented a first evi- 
dence of D° mixing in decays of charmed 
mesons to CP eigenstates [12j and the most 
precise measurement of the mass difference 
in the neutral charmed meson system [lq ]. 
The 68% confidence regions of mixing pa- 
rameters x and y arising from the two mea- 
surements are presented in Fig. [9l 
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01 
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FIG. 9: 68% confidence level regions of x and 
y arising from the lifetime measurements in 
D° —¥ K + K~, tv + tt~ decays and time depen- 
dent Dalitz analysis of D° — > Ksir + iv~ . 



The charm subgrou p o f the Heavy Fla- 
vor Averaging Group [2CJ has at the time 
of the conference presented world averages 
of measurements in the field of charm mix- 
ing and CPV. By summation of the like- 
lihood curves depending on various observ- 
ablcs (and thus accounting for non-Gaussian 
distribution of some experimental uncertain- 
ties) the n a (n = 1 — 5) 2-dimensional con- 
tours in (x, y) plane are presented in Fig. 1101 
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0.03 
0.02 
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0.04 -0.03 -0.02 -0.01 0.01 0.02 0.03 0.04 

x 

FIG. 10: Contours of n a 2-dimensional allowed 
regions for x and y from the world average of 
measurements assuming no CPV . 

The shape of the contours is a non-trivial 
result of summing likelihoods of specific 
shapes. Especially the constraints arising 
from the D° — > K + it~ decays, when ro- 
tated in (x,y) plane due to the uncertainty 
in knowledge of the strong phase difference 5, 
exhibit an almost circular region of increased 
values around the (0, 0) point. 



* = (o.87 ± °;>„ 

21 

y = (0.66± Q20 )% 



(19) 



The no-mixing point (x, y) — (0,0) is ex- 
cluded at more than 5 standard deviations. 

Due to a larger number of free parame- 
ters the method of averaging experimental 
results which allow for a possibility of CPV 
is a x 2 minimization. Nevertheless the re- 
sults for the mixing parameters are almost 
unchanged, and CPV parameters are consis- 
tent with no violation of the CP symmetry: 



x = (0.84 ± 



0.32 N 



0.34' 
y = (0.69 ±0.21)% 

M=o.88 ± °- 23 
0.20 



(-0.09 ± 



0.17, 
0.19' 



rad 



(20) 



The world averages are dominated by re- 
sults from the existing B-factories, with a sig- 
nificant contribution from Cleo-c and Teva- 
tron 2l|. The results presented in this pa- 
per were obtained by the Belle collabora- 
tion using around one half of the expected 
full data set. The experimental errors are 
mainly dominated by statistical uncertain- 
ties and will thus improve in the next year. 
However, considering the uncertainty of the 
SM predictions as well as the small CPV 
expected within, it is unlikely that the full 
range of high scientific interest in the mea- 
surements of the charm mixing will be ful- 
filled by the end of the data taking of Belle 
and BaBar. It is thus instructive to make an 
attempt of predicting the accuracy that may 
become available at some of the future exper- 
iments. In Table |TT] the expected one stan- 
dard deviation errors on the key parameters 
are given for an average of measurements to 
be performed at the proposed Supcr-B fac- 
tory. The values are estimated by scaling 
the current Belle statistical sensitivity and 
an educated guess on possible improvements 
of systematic uncertainties. Two values of 
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TABLE II: Expected one standard deviation er- 
rors on the measurements of charm mixing and 
CPV parameters, to be performed at the Super- 
B factory. 



Parameter 


C = 5 ab^ 1 


C = 50 ab^ 1 


a(x) 


0.1% 


0.07% 


°(y) 


0.1% 


0.07% 


a{A M ) 


0.1 


0.07 


o{4>) 


0.1 rad 


0.07 rad 



a precise determination of x and search for 
the CP violation in the charm sector. For 
it is in this field of measurements where we 
can expect an answer to the above question. 
Observation of CPV effects at the existing 
facilities would be a sign of NP. A proposed 
Super-B factory would enable searches for 
the CP asymmetries to the 10 -4 level, a 
range covered by the current SM expecta- 
tions and thus an interesting area of search 
for NP to appear. 



expected integrated luminosity are consid- 
ered; the lower one represents a modest data 
sample that could be collected at the Super- 
B factory while the higher one is an ultimate 
goal. The a's shown are of course to be taken 
with a grain of salt since the systematic er- 
rors which are difficult to estimate are impor- 
tant if not the dominating part of the total 
error. Nevertheless they are worth present- 
ing, if for nothing else then for an easy-to- 
remember pattern. 

The expected sensitivities can be com- 
pared to expectations from the LHCb exper- 
iment scheduled to start data taking in the 
coming year. With C — 10 fb" 1 (expected 
to be collected in around 5 years of run- 
ning at the nominal luminosity) one hopes 
for a(x' 2 ) « 6 x 10Z\ <r(y') « 0.9 x 10~ 3 and 
&{ycp) ~ 0.05% |22j . These estimates can 
be roughly placed in between C = 5 ab _1 and 
C = 50 ab _1 expectations from the Super-B 
factory. 

Considering the results presented in this 
and other charm mixing related papers sub- 
mitted to the conference, it is fair to say that 
the year 2007 (31 years after the D° meson 
discovery) was the year of experimental con- 
firmation of the D° mixing (to be compared 
to the time span of 6 years between simi- 
lar observations in K° system, 4 years in 
system and 14 years in the case of B® 
mesons). At the moment we are facing a 
somewhat rare situation of experimental ev- 
idence without an accurate theoretical guid- 
ance on whether the phenomenon is entirely 
due to the SM physics or not. The largest 
interest determining the work ahead lies in 
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